Ti-enhanced kinetics of hydrogen absorption and desorption on NaAlH4 surfaces 
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We report a first-principles study of tfie energetics of iiydrogen absorption and desorption (i.e. 
H-vacancy formation) on pure and Ti-doped sodium alanate (NaAlH4) surfaces. We find tiiat 
the Ti atom facilitates the dissociation of H2 molecules as well as the adsorption of H atoms. In 
addition, the dopant makes it energetically more favorable to creat H vacancies by saturating Al 
dangling bonds. Interestingly, our results show that the Ti dopant brings close in energy all the 
steps presumably involved in the absorption and desorption of hydrogen, thus facilitating both and 
enhancing the reaction kinetics of the alanates. We also discuss the possibility of using other light 
transition metals (Sc, V, and Cr) as dopants. 

PACS numbers: 81.05.Zx, 81.05.Je, 61.12.-q, 63.20.Dj 
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In recent years sodium alanate (NaAlH4) has become 
one of the most promising systems to achieve safe and 
inexpensive storage of hydrogen aboard vehicles. The ki- 
netics of the reversible reaction by which pure NaAlH4 
releases and stores hydrogen is relatively slow. However, 
it was discovered that a few percent of Ti doping in- 
creases the reaction rates dramatically [l|, bringing the 
system close to what is required for practical applica- 
tions. In spite of extensive investigations, the mecha- 
nism by which Ti enhances the reaction kinetics remains 
largely unknown Q,y,|j]. As a matter of fact, only re- 
cently have we had convincing experimental 5] and the- 
oretical [a, [3 evidence that the Ti dopants are located 
on the surface, as opposed to the bulk, of the system. 
First-principles simulations by several groups also pre- 
dict that Ti dopants cause the breaking of multiple Al-H 
bonds H, Hi and that bulk Ti dopants |3| and Na va- 
cancies [13 facilitate the creation of H vacancies. Such 
studies suggest mechanisms by which Ti can improve the 
H desorption kinetics of the alantes, but are far from be- 
ing complete or conclusive (specially those restricted to 
the bulk of the system). 

In this work we present a first-principles study of the 
energetics of hydrogen absorption and desorption (i.e. 
vacancy formation) on pure and Ti-doped NaAlIl4 su- 
faces. We find that the Ti dopant facilitates both the 
absorption and the desorption, thus enhancing the cy- 
cling kinetics of hydrogen. 

The calculations were performed within the general- 
ized gradient approximation [ij to density functional 
theory [12 as implemented in the code SIESTA [l^. 
We used a localized basis set including double-^ and po- 
larization orbitals, and TrouUier-Martins pseudopoten- 
tials jl4j • We tested the convergence of our calculations 
with respect to the k-point and real space meshes: we 
typically considered a 4500 A'^ simulation supercell con- 
taining about 162 atoms, for which we used a 2x2x1 
k-point grid and a 150 Ry cutoff for the real space mesh. 
The accuracy of our SIESTA resutls had been previously 
tested against more accurate ah initio calculations using 
a plane- wave basis set Q- 

Tetragonal NaAlH4 presents two natural terminations 



related to the (100) and (001) directions, respectively. 
The phenomenology associated with the Ti-doping of 
these two surfaces is known to be essentially the same V] ; 
thus, for this work we considered the 001-surface, which 
can be treated using a relatively small supercell. Our 
slab-type supercell can be regarded as composed of 
3x3x3 NaAlH4 groups, i.e., 162 atoms. We imposed 
that the atoms in the deepest layer be fixed at the 
bulk atomic positions, and allowed about 13 A of empty 
space between slabs. The bulk structural parameters 
(a = b =5.01 A and c =11.12 A) were obtained from first- 
principles 13- Structural relaxations were considered 
to be converged for residual force components smaller 
than 0.02 eV/A. We calculated absorption and vacancy- 
formation energies as the difference between the energies 
of products and reactants; thus, a negative value signals 
an energetically favorable reaction. To study, for exam- 
ple, the deposition of H atoms, we considered several 
starting configurations, performed regular structural re- 
laxations, and took the lowest-energy result. Ideally one 
would use, for example, simulated annealing techniques 
to identify the true lowest-energy structures, but such an 
approach is too computationally costly. 

The question of whether the Ti dopant substitutes for 
Na or Al is somehow controversial, as it critaically de- 
pends on the choice of energies that are used as a ref- 
erence to compute the relative stability of the various 
doping models [ifl ■ The Ti doping of the alanates is usu- 
ally done by ball-milling, which is an out-of-equilibrium 
process that involves individual Ti, Na, and Al atoms, as 
well as H2 molecules, rather than the corresponding crys- 
talline phases. Hence, we prefer using atomic energies as 
reference, which leads to the conclusion that Ti substi- 
tutes for Na 13, Hi- However, if the energies of crystalline 
Ti, NajEind Al are used, the result is that Ti substitutes 
for Al 6]. Being aware of this issue, in this work we have 
considered these two Ti-doping models, i.e. Ti substitut- 
ing for Na (Ti@Na) and Ti substituting for Al (Ti@Al), 
and found that both lead to the same qualitative results. 
Note also that our present discussion is based on how 
the pure and doped surfaces differ regarding the energet- 
ics of hydrogen absorption and vacancy formation. Our 



TABLE I: Absorption energies (in eV) of molecular H2 and 
various nH cases (where the H atoms are deposited on the 
surface one at a time) of pure and Ti-doped NaAlH4 surfaces. 
Two doping models are considered, in which Ti substitutes 
for Na (Ti@Na) and Al (Ti@Al) , respectively. We use as ref- 
erence the energy of the surfaces with no extra hydrogen, and 
the energy of the H2 molecule. See text for the explanation 
of the Ti®Na-|-H2 results. 



Surf, type 



+H2 



-H 



+2H +3H +4H 



pure 


-0.07 


1.81 






Ti@Na 


-0.03/-0.23 


-0.49 


-0.23 


-0.67 


Ti®Al 


-0.41 


-0.61 


-0.41 


-0.60 



-0.38 



conclusions are thus independent of the chosen energy 
reference, which is the same for both pure and doped 
surfaces. 

Hydrogen absorption.- We first studied the possibil- 
ity that the alanate surface, pure and Ti-doped, might 
absorb extra hydrogens. The underlying assumption is 
that a favorable absorption energetics will facilitate the 
NaAlH4-formation reaction, which occurs under H2 pres- 
sure. We considered the direct absorption of an H2 
molecule as well as the absorption of individual H atoms 
that we deposited on the surface one at a time. Tabled 
summarizes our results. 

We found that the surface is able to bind an H2 
molecule in all the cases considered. The molecule disso- 
ciates in the Ti@Al case, forming an AlHg group, which 
renders a significant energy gain of 0.41 eV. In the case 
of the pure surface, the H2 molecule does not dissociate 
and the binding energy is significantly smaller (0.07 eV). 
Finally, the Ti@Na case turns out to be more compli- 
cated: When the H2 molecule is deposited on the relaxed 
Ti@Na surface, it does not see the Ti atom and binds 
weakly without dissociating (the calculated binding en- 
ergy is 0.03 eV). However, if the H2 molecule is deposited 
on the unrelaxed Ti@Na surface, it dissociates and binds 
with an energy of 0.23 eV. It is useful to think of this 
absorption as the two-step process depicted in Fig.^ In 
the first step the molecule breaks and the H atoms bind 
to the Ti on the alanate surface. The corresponding en- 
ergy gain is about 1.4 eV, which is very close to what 
is obtained for H2 breaking and absorption on carbon 
nanostructures coated with Ti [13. In the second step 
the whole surface relaxes and the Ti dopant adopts a fi- 
nal configuration with 8 neighboring H atoms (see botom 
inset of Fig. ^ . In conclusion, our results for both Ti@ Al 
and Ti@Na clearly show that the Ti dopand facilitates 
H2 dissociation and absorption. 

Next we studied the deposition of individual H atoms 
one at a time. We found that the pure surface is unable 
to absorb an individual H atom (note the positive absorp- 
tion energy of 1.8 eV in Tabled; in fact, the deposited H 
destabilizes an AIH4 group and forms an H2 molecule. In 
contrast, the doped surfaces are able to absorb the first 
H atom with significant energy gain (about 0.5 eV); the 
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FIG. 1: (color online) Energy versus structural relaxation step 
for the H2 dissociation and absorption by the TiQNa surface. 
The starting point is the unrelaxed Ti@Na surface with an 
H2 molecule close to it, which is taken as the zero of energy. 
In the first part of the relaxation, only the H2 molecule is 
allowed to move; then, all the atomic positions are relaxed. 
Given for reference is the result for the relaxed Ti@Na plus 
an H2 located at infinity. The insets show the various steps 
of the absorption process. 



H atom locates in the vicinity of the Ti dopant, at about 
1.9 A, forming a TiH5 group in the case of Ti@Al. More- 
over, further hydrogens can be absorbed by the doped 
surfaces. For Ti@Na, the -I-3H situation seems to be par- 
ticularly stable; for +4H an II2 molecule forms and stays 
bound to the surface, but not properly absorbed by it. 
For Ti@Al, we found that a H2 molecule forms already 
in the -t-3H case. Thus, our results indicate that a Ti 
dopant can take a maximum of 3 or 4 extra hydrogens. 

In order to enhance the NaAlIl4-formation kinetics, 
the absorbed hydrogens should be mobile on the sys- 
tem surface. To check this, we performed computer ex- 
periments as the following one: we start from the re- 
laxed Ti@Na-|-2H system, create a H vacancy and thus 
a defective AIH3 group, relax the system (H migration 
should occur to recover an AIH4 group), and compare 
the obtained structure and energy with the assumed 
lowest-energy solution Ti@Na-|-H. We observed two gen- 
eral trends. (1) Large structural distortions and atom 
migrations occur. Most relaxed structures differ signif- 
icantly from the lowest-energy solution and are slightly 
above in energy (by about 0.1 eV, typically). (2) Hydro- 
gen migration occurs more easily in H-rich surfaces, i.e., 
in the -I-3H and -I-4H cases of Table U The difficulties 
to reach the reference, lowest-energy structure reflect the 
fact that we are dealing with a complicated multi-minima 
energy landscape. On the other hand, the significant 
atomic mobility observed suggest that, at room temper- 
ature, the H atoms initially bound to the Ti dopant will 



TABLE II: Desorption energies (in eV) of hydrogen and other 
atom groups for pure and two Ti-doped NaAlIl4 surfaces. We 
denote by +H" the surface with one H vacancy, etc. We use 
as reference the energy of the original pure and Ti-doped sur- 
faces, the energy of the H2 molecule, and the atomic energies 
of Na and Al. 
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Surf, type 


-HH" 


-f-Na" 


-H(NaH)" 


-KAIH4)" 


1^ 


pure 


1.98 


7.96 


5.61 


7.51 




Ti@Na 


0.41 


4.89 


5.41 


5.70 


t§ 


Ti@Al 


0.30 


5.20 


4.50 


5.65 





indeed migrate and take part in the NaAlH4-forniation 
reaction, thus increasing its rate. 

Hydrogen vacancy formation.- Next we studied the 
ability of the Ti dopant to facilitate the formation of 
hydrogen vacancies and, further, destabilize the surface. 
We calculated the formation energies of various types of 
vacancies: H, Na, and atom groups like NaH and AIH4. 
(We denote the H vacancy by H", etc.) Table Hll summa- 
rizes our results. 

We first note that all the vacancy-formation energies 
are positive; thus, as expected, wc find that the surface, 
either pure or doped, is stable and thermal activation en- 
ergy is required to decompose it. The second observation 
is that the pure surface is significantly more stable than 
the Ti-doped surfaces. For example, the formation en- 
ergy of H vacancies is about 2.0 eV for the pure surface, 
while it is below 0.5 eV when the Ti dopant is present. 
The same trend is observed in the formation energy of 
the other studied defects, which clearly indicates that the 
Ti-doped surface will decompose at temperatures signif- 
icantly lower than those required in the case of a pure 
system. 

The Ti dopant is able to reduce the defect formation 
energies thanks to its chemical versatility. In the case of 
charged defects, the dopant adquires extra charge (e.g. 
in the Ti@Na-|-(AlH4)'' case) or enables effective passi- 
vation of dangling bonds. The latter is best illustrated 
by the Ti@Al-f ff' case depicted in Fig. |2 The relaxed 
configuration presents the Ti and a neighboring Al atom 
within bonding distance (2.62 A) and sharing two hydro- 
gen atoms. This relaxation results in a significant en- 
ergy reduction and thus facilitate the H" formation. In 
the case of neutral defects (e.g. (NaH)") the Ti dopant 
does not seem to play such a relevant role and, thus, 
the pure and doped surfaces present relatively similar 
vacancy-formation energies. 

Note that we considered explicitly the case of Na- 
vacancy mediated H desorption, which has recently been 
claimed to be specially important on the basis of bulk 
NaAlIl4 calculations [l^- Indeed, we found that Na- 
defective pure and Ti@Al surfaces gain energy by releas- 
ing a hydrogen atom (about 2.3 and 0.7 eV, respectively), 
while there is an energy cost of 0.5 eV in the Ti@Na case. 
On the other hand, the occurrence of the Na vacancy is 
much more likely when the Ti dopant is present, for both 




50 100 

Minimization Iteration (steps) 



150 



FIG. 2: (color online) Energy versus structural relaxation step 
for the relaxation of a Ti@Al surface in which we create a hy- 
drogen vacancy neighboring the Ti atom (we denote this case 
by Ti@AH-H"). The formation of an Al-Ti bond in the final 
configuration results in a reduced H" formation energy. The 
insets show the initial and final configurations of the atoms 
mainly involved in the structural relaxation. 



Ti@Al and Ti@Na, which renders the net (NaH)" forma- 
tion energies in Table ^ As already mentioned, since 
(NaH)" is not a charged defect, there is not a very big dif- 
ference between the results for pure and doped surfaces. 
Our calculations indicate that the Ti dopand favors the 
Na"-mediated mechanism for H desorption slightly in the 
Ti@Na case and more significantly in the Ti@Al case. 

Our calculations render a conclusion that may seem 
paradoxical at first, namely, that the Ti-dopant fa- 
vors both the absorption and desorption of hydrogen on 
NaAlH4 surfaces, thus enhancing the kinetics of the H- 
charge and discharge reactions of the system. The expla- 
nation to this paradox is that, as a result of its chemical 
versatility, the Ti dopant brings close in energy all the 
steps presumably involved in these reactions, with the 
corresponding reduction in activation energies and, thus, 
increased reaction rates. For this picture to hold, we need 
to assume that the Ti dopants remain on the surface of 
the system during the reactions, which seems plausible 
given the large ionic mobility involved in these processes 
and the fact that, as shown in Ref. |3], it is energeti- 
cally favorable for the dopant to lie on the surface of the 
system. 

Finally, we studied the possibility of using, instead of 
Ti, other light transition metals as dopants. More pre- 
cisely, we tested the performance of Sc, V, and Cr by 
calculating the energies of absorption and desorption of 
a single H atom by a Na-substituted surface. Regarding 
absorption, only the Sc dopant favors it, with a negative 
absorption energy of —0.27 eV; we get 0.26 and 0.20 eV 
for V and Cr, respectively. We also found that Sc causes 



the H2 molecule to dissociate provided H2 is deposited 
on the unrelaxed Sc@Na structure. As for the formation 
of vacancies, the lowest energy is obtained for Cr with 
0.30 eV; V gives 0.75 eV and Sc renders the most sta- 
ble surface, with a formation energy of 1.13 eV. Hence, a 
comparison with the pure surface results in Tables Q and 
^indicates that Sc, V, and Cr will facilitate both hydro- 
gen absorption and desorption. At the same time, our re- 
sults clearly show that Ti performs better in all aspects. 



thus being the dopant of choice among light transition 
metal atoms. 
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